arXiv:1502.05687v3 [cond-mat.mtrl-sci] 23 Nov 2016 


Characterization of spin relaxation anisotropy in Co using spin pumping 

Y. Li, W. Cao and W. E. Bailey 

Materials Science & Engineering, Department of Applied Physics and Applied Mathematics, 
Columbia University, New York NY 10027, USA 
(Dated: November 24, 2016) 

Ferromagnets are believed to exhibit strongly anisotropic spin relaxation, with relaxation lengths 
for spin longitudinal to magnetization significantly longer than those for spin transverse to mag¬ 
netization. Here we characterize the anisotropy of spin relaxation in Co using the spin pumping 
contribution to Gilbert damping in noncollinearly magnetized Pyi_a;Cua;/Cu/Co trilayer structures. 
The static magnetization angle between Pyi_a,Cua; and Co, adjusted under field bias perpendicular 
to film planes, controls the projections of longitudinal and transverse spin current pumped from 
Pyi-ajCux into Co. We find nearly isotropic absorption of pure spin current in Co using this tech¬ 
nique; fits to a diffusive transport model yield the longitudinal spin relaxation length < 2 nm in Co. 
The longitudinal spin relaxation lengths found are an order of magnitude smaller than those deter¬ 
mined by current-perpendicular-to-planes giant magnetoresistance measurements, but comparable 
with transverse spin relaxation lengths in Co determined by spin pumping. 


A key question for spin electronics concerns the relax¬ 
ation mechanisms for spin current injected into a vari¬ 
ety of materials. Spin relaxation in ferromagnets (Fs), 
central for spin momentum transfer, is special because of 
the anisotropy axis presented by the spontaneous magne¬ 
tization M |ll-[lol|. Longitudinal spin relaxation[i|, with 
spin polarization cr parallel (antiparallel) to M, causes 
spin accumulation to decrease exponentially with dis¬ 
tance over a scale greater than the electronic mean free 
pathQ. Transverse spin relaxation, with a orthogonal 
to M, is governed by the dephasing process of spin-up 
and spin-down eigenmodes due to their different Fermi 
wavevectors, leading to oscillation and decay of spin ac¬ 
cumulation on a scale shorter than the electronic mean 
free pathjl, [I^. 

The characteristic length scales for the two different 
spin relaxation processes in ferromagnets, for longi¬ 
tudinal and for transverse spin relaxation, have been 
evaluated largely using two separate experimental tech¬ 
niques: magnetotransport [Ilj for A^ and ferromagnetic 
resonance (FMR)3,[i] for A^. These two measurements 
characterize charra-accompanied and chargeless spin cur¬ 
rent, respectively[l|,lj|. Estimates of A^ come from the F 
layer thickness dependence of current-perpendicular-to- 
planes giant magnetoresistance (CPP-GMR) ; ex¬ 

tracted values of A^j. range from 5 nm for Ni 79 Fe 2 i up to 
40 nm for Co at room temperature. FMR measurements 
of spin pumping, for collinearly magnetized Fi/Cu/F 2 
structures, show much shorter penetration depths (Ac) 
to fully absorb transverse spin current [H El • Co has 
the most anisotropic spin relaxation according to these 
separate measurements, with Af;./A^ ~16 taking A^ ^ 
2Ac = 2.4 nmEl,E3- 

In this manuscript, we demonstrate that the longitu¬ 
dinal spin relaxation length, in addition to the trans¬ 
verse spin relaxation length [l5|, can also be character¬ 
ized using a spin pumping measurement, enabling a mea¬ 
surement of the anisotropy of spin relaxation in a given 


ferromagnetic layer. We present FMR measurements 
of the spin pumping contribution to Gilbert damping 
in noncollinearly magnetized Pyi-^Cux/Cu/Co multi¬ 
layers (Py=Ni 7 gFe 2 i). Using Pyi_a;Cua; alloys, which 
have adjustably smaller saturation magnetization Ms 
than Co, we can change the magnetization alignment of 
Pyi_a;Cua; and Co from collinear for in-plane FMR to 
near-orthogonal for perpendicular FMR. As the angle 0 m 
between Pyi_a;Cua; and Co magnetization tends towards 
7r/2, one component of injected spin from Pyi_ 3 ;Cua; 
tends towards the longitudinal direction (Fig. 1), allow¬ 
ing us to probe anisotropy in spin relaxation through the 
linewidth of the Pyi-^Cux layer 18, Ei- We find, sur¬ 
prisingly, that spin relaxation, as measured through the 
spin pumping contribution to Gilbert damping, is mostly 
isotropic. In our Go hlms we estimate A^ < 2 nm for all 
different Pyi_ 2 ;Gua;/Gu/Go samples, which is compara¬ 
ble to its transverse counterpart(^2.4 nm) but inconsis¬ 
tent with the much longer (^40 n m) lengths reported 
from room-temperature GPP-GMR[iiI IT3| . 

Three types of thin-film heterostructures were pre- 



FIG. 1. Left: Noncollinear magnetization alignment of the 
F 1 /N/F 2 trilayer at the FMR condition for Fi. Right: mi 
is driven into precession, pumping spin current into m 2 , with 
spin components both longitudinal (cl) and transverse (ct) 
to the m 2 magnetization. 

pared by UHV sputtering and characterized by FMR. 
Pseudo-spin-valve-type Pyi_a;Cua;(t)/Cu(5 nm)/Co(5 
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nm) trilayers were used to characterize the anisotropy 
of spin-current absorption in Co. Their response was 
compared with two types of Pyi_ 2 ;Cu 2 ;(t) alloy control 
samples. Bilayers of Pyi_a;Cua;(5 nm*)/Cu(5 nm) and 
trilayers of Pyi_a;Cua:(5 nm*)/Cu(5 nm)/Pt(3 nm) were 
used to characterize the background damping of the alloy 
and the spin mixing conductance of the alloy/Cu inter¬ 
face, respectively. Co(5 nm)/Cu(5 nm)/Pt(3 nm) is also 
deposited. For the alloy Cu contents x = 0 to 0.4 were 
prepared in each case, using confocal sputtering from Py 
and Cu targetsjl^; thicker (10 nm*) alloy layers were 
used for x = 0.4. All layers were deposited on Si/Si 02 
substrates, seeded by Ta(5 nm)/Cu(5 nm) and capped 
by Ta(2 nm). See Ref. [2lj for details on preparation. 

Room temperature, variable frequency (3-26 GHz), 
swept-field FMR measurements were used to character¬ 
ize the samples, with instrumentation as described in 
(^ . In order to characterize FMR relaxation of the 
Pyi-xCua, layer under noncollinear magnetization align¬ 
ment with Co, two types of measurements were car¬ 
ried out. First, we compare the frequency-dependent 
linewidths of Pyi_a;Cux and Pyi_a;Cua;/Cu/Co sam¬ 
ples in both in-plane (parallel-conditio n, p c) and per¬ 
pendicular (normal-condition, nc) FMR[23j, for a series 
of four measurements at a given alloy content x\ see 
Figs. Et, E and m Here we expect the Co magneti¬ 


zation of trilayer samples to vary from fully perpendic¬ 
ular to the film plane at high biasing field (high 
Loj^'K) to nearly parallel to the film plane at low iJe 
(low while the Pyi_a;Cua; magnetization is al¬ 

ways perpendicular to the him plane. Second, we com¬ 
pare the polar angle-dependent linewidths of Pyo. 8 Cuo .2 
and PyQ. 8 Cuo, 2 /Cu/Co samples at a Hxed frequency of 
a;/27r = 10 GHz; see Fig. [S] Here we expect the mis¬ 
alignment angle 0 m to change from zero to maximum as 
we rotate the biasing Held from in-plane (pc) to out-of- 
plane (nc). 

Theoretical models for the spin pumping contribution 
to damping under noncollinear magnetization alignment 
of symmetric F \ /N/Fi structures were developed in Refs. 
[3 [13. We have extended these models to consider 
asymmetric F 1 /N/F 2 structures where Fi=Pyi_ 2 ,Gua; 
and F 2 =Go in our samples. In the spin valve structure 
the spin-pumping damping enhancement Aogp of Fi is 
caused by the dissipation of spin current in F 2 . If Fi and 
F 2 are misaligned by an angle 0 m, where 0 m = mi • m 2 
(Fig. [T|), during small-angle precession of Fi, the polar¬ 
ization of spin current pumped into F 2 will oscillate from 
fully transverse to maximally longitudinal. The instan¬ 
taneous spin-pumping damping will then oscillate from 
asp( 0 °) = A On X / ( 5 i'*' + 52 ^), as given in the standard 
collinear casejl^, to a minimum value given bydl: 


Aasp(0M) = Atto 


^2 (A sin^ 0M — H sin 0 m cos 0 m ) + gV'iC cos^ 0 m — H sin 0 m cos 0 m ) 


AC-B^ 


( 1 ) 


Here and g* {i = 1, 2) are the effective trans¬ 
verse and longitudinal spin conductances, respectively; 
Aoq = 7?i5i'*'/(47rMsfi?) is the damping enhancement 
with effective spin mixing conductance of 
in Eq. [T] A(0 m) = 5 * sin^ 0 m + cos^ 0 m -b ffl'*', 
B{0m) = (gl'‘' - 5i ) sin 0m cos 0M and C(0m) = 

gl cos^ 0M + 1 '*' sin^ 0 m + 52 ■ We take the arithmetic 
mean of the two extreme cases as the effective damping 
enhancement, as found to be valid in Ref. El- See the 
Supplemental Materials for details. 

To maximize the spin pumping anisotropy at Hnite 
0M, we use Go (5 nm) for F 2 , where the dimension is 
chosen to be signiHcantly thicker than the transverse 
spin penetration depth, Ac = 1-2 nm 15l|. and thinner 
than the reported longitudinal relaxation length A^, 
^ 38 nm[l3, resulting in a large expected asymmetry 
in spin relaxation. In the analysis of relaxation in non¬ 
collinear ly magnetized structures, we take spin mixing 
conductances g/^^ as parameter inputs, determined from 
the measurements on the Ft control structures, and 
take the longitudinal spin relaxation length A^j. as a Ht 
parameter. 


Fig. E] summarizes the results of Hxed-angle nc-and 
pc-FMR measurements for the three sample series. In 
Fig. 2(a) we plot resonance Helds noHres as a function of 
frequency for single layers and trilayers in nc-FMR. The 
good agreement in the ^o-Rres of Pyi-xCua, measured in 
single layers and trilayers demonstrates that Pyi_:jGua; 
properties are reproducible in deposition. In Fig. 2(b) 
the effective magnetizations g-oMeg, extracted from Hts 
to the linear Kittel equation w /7 = /iq (Fires ~ are 

plotted as a function of x. The data show Slater-Pauling 
dilution of magnetic moment in the Pyi-a^Gua, layer 
with increasing Gu content 

In Fig. Elc) we plot full-width half-maximum linewidth 
goAHi /2 as a function of uj/27t at a;=0.2. Gilbert-type 
damping, /j,oAiFi /2 = MoAiFo + with negligible 

inhomogeneous broadening ^qAIFo, is observed for both 
pc- and nc-FMR in the single layer and for pc-FMR in 
the trilayer. The linewidths agree closely for pc- and 
nc-FMR in the single layer, showing a ne glig ible role 
for two-magnon scattering in the linewidth|26l|. In the 
trilayer, nc and pc linewidths agree well for frequencies 
above 10 GHz. These observations hold for samples 
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FIG. 2. (a) Perpendicular (nc-FMR) resonance field fioHiea 
for Pyi_a;Cua; single layers and Pyi_a;Cua;/Cu/Co trilayers, 
a; = 0 - 0.4, as a function of frequency a;/27r. (b) Effective 
magnetization fioMes extracted from (a) as a function of x. 
(c) Resonance linewidtlis fioAHi /2 of the Pyo.gCuo.z single 
layer and trilayer as a function of frequency lijI2h. The spin 
pumping enhancement is clearly visible in the increased slope 
(a) of the trilayer data; the low-frequency deviation is dis¬ 
cussed in Fig. 3. (d) Effective spin mixing conductances 
of Pyi_a;Cua;/Cu/Co, Pyi-xCua;/Cu/Pt and Co/Cu/Pt. 

with all Cu content 0 > x > 0.4; the deviations at 
low frequency are discussed in Fig. [31 The effective 
spin mixing conductances gjg of trilayer samples are 
extracted from Aosp = 7^ffJff/(47rMstM), shown above 
where Aosp is the difference in a between trilayers and 
single layers. In Fig. 2(d) we show the extracted 
for Pyi_a;Cua;/Cu/Co and Pyi-xCux/Cu/Pt structures 
as a function of x. We also plot gjg of Co/Cu/Pt 
for reference. The lower level of gjg ^7 nm“^ for 
Pvi -xCux/Cu/Co, compared with ^15 nm“^ measured 
in Ref. [l5[, is likely to be from a more resistive Cu layer, 
which adds an additional resistance of {2e^/h)tcu/<xcu 
to the inverse of total spin mixing conductance where 
CTCu is the Cu conductivity. Using these gig values, 
we extract the effective spin mixing conductance of 
Pyi-xCux/Cu and Co/Cu interfaces, shown in the 
Supplemental Materials [2^. These parameters will 
be used to determine the longitudinal spin relaxation 
lengths from the spin pumping data in Figs. 0]and|5] 

In the measurements presented in Fig. |2lc), the nc- 
FMR linewidths are measured at applied fields below the 
saturation field for Co, fioMes = 1-4 T. The saturation 
field corresponds to a nc-FMR resonance frequency for 
Pyo, 8 Cuo .2 of 25 GHz, as shown in Fig. UKa). With the 
resultant noncollinear magnetization alignment in the 
trilayer, we expect to see spin-pumping damping Aogp 
for Pyo. 8 Cuo ,2 reduced in nc-FMR compared with the 
values in pc-FMR. Instead, we find that the linewidths 
of the trilayer measured in pc- and nc-FMR agree closely 
when a;/27r > 10 GHz. Furthermore, in nc-FMR there is 


an additional broadening from 2-10 GHz in the trilayers 
which is not predicted by the model. 

In order to determine whether the low-frequency 



FIG. 3. pc- and nc-FMR linewidths for single (Py) and 
trilayer (Py/Cu/Co) structures, introducing MgO interlayers 
to suppress spin pumping. Dashed lines are linear fits to pc- 
FMR linewidths. Solid curves assume (magnetostatic) inter¬ 
layer coupling of 10 mT acting on Py and reproduce the low- 
frequency upturn in linewidth, seen to be present equally with 
and without MgO. Inset: enhancements of nc-FMR linewidth 
over pc-FMR linewidth for the three samples. 

broadening is related to spin pumping, we have also 
measured pc- and nc-FMR linewidths of Py(5 nm)/Cu(5 
nm)/MgO(2 nm) and Py(5 nm)/Gu(3 nm)/MgO(2 
nm)/Co(5 nm) structures, deposited with the same 
seed and capping layers. MgO interlayers are known to 
suppress spin pumping (^. Introducing MgO between 
Py and Co, we show in Fig. [3] that the pc linewidths of 
Py in trilayer Py/Cu/Co (blue crosses) are restored to 
those of single-layer Py/Cu/MgO (overlapping green and 
red crosses), demonstrating suppression of spin pumping 
between Py and Co. However, we see a very similar 
upturn in low-frequency (< 10 GHz) Py linewidth in 
nc-FMR (red circles), similar to that shown in Fig. 
He). We attribute this low-frequency behavior to an 
interlayer coupling which cants the magnetization of 
Py a few degrees off the film normal when Co is not 
fully saturated along the film normal (i.e. Hq < Mes). 
The solid curves in Fig. 3 assume a coupling field of 
10 mT on Py, parallel to the local Co magnetization, 
which reproduce the linewidth broadening of nc-FMR. 
The peak-like features around 3 GHz show the maximal 
Gilbert damping enhancement when the Py magnetiza¬ 
tion is canted, as demonstrated in Fig. 5 inset. 

Fig. 4 shows the central result of the paper. We 
compare the spin-pumping linewidth enhancements, 
/io(Ai/*() 2 ~^^r/?*'*)’ between pc- and nc-FMR (crosses 

and circles) in Fig. 4(a-d). Here and AFIj ^"2 

are the linewidths of Pyi-xCux in Pyi-^Cux/Cu single 
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layers and Pyi_a;Cua;/Cu/Co trilayers, respectively. The 
spin pumping linewidths are quite linear as a function of 
frequency for the pc-FMR data, as expected. However, 
above 10 GHz (shaded regions), they are also quite 
linear in nc-FMR, which is not expected. Collinear 
and noncollinear spin pumping linewidths agree closely. 
This behavior is in contrast to the predicted behavior 
using A^j. = 38 nm for Co, measured by CPP- GMRini, 
and calculated in dashed curves according to the theory 
in the Supplemental Materials. From the evident 
agreement between pc- and nc-linewidths above 10 GHz, 
for all Cu content x, we find no evidence for anisotropy 
in spin relaxation in our Co films. Best fits to the data 
yield longitudinal spin relaxation lengths < 2 nm 
in each of the four cases, approximately equal to the 
previously measured transverse length A^ = 2.4 nm[l5|. 

Our model has assumed single-domain (macrospin) 


X pc-FMR • nc-FMR 



FIG. 4. Spin pumping contribution to linewidth in pc- and 
nc-FMR. (a-d) Linewidth enhancement of Pyi_a,Cua; between 
single layers and trilayers in pc- and nc-FMR, x = 0.1-0.4. 
Solid lines are linear fits to the pc data (crosses); dashed 
curves are predicted from Eq. (1) using A^=38 nm. The 
shadows at a;/27r < 10 GHz denote where the low-frequency 
linewidth broadenings are significant. 

behavior in both Co and Pyi_a;Cua; layers. For 
Pyi_a;Cua; under field bias well in excess of M^, the 
magnetization is well saturated, but for the Co layer, 
with higher Mg, nonuniform magnetization is possible. 
For greater control over the Co domain state, we have 
also carried out angle-dependent, fixed-frequency FMR 
measurements on Pyo. 8 Cuo .2 and Pyo, 8 Cuo, 2 /Cu/Co. 
Here the Co layer can be saturated more easily be¬ 
cause the biasing field is canted away from the normal 
condition. The frequency is set to 10 GHz, where the 
low-frequency linewidth broadening of Pyo. 8 Cuo ,2 is 
insignificant (Figs. 3 and 4). As the field angle 0h 
goes from 90° to 0° (pc to nc), the angle between the 
magnetizations of Pyo. 8 Cuo .2 and Co changes from zero 


to maximum noncollinearity (^ 50°) and Aogp would 
be expected to decrease significantly where the spin 
relaxation length in Co is markedly anisotropic. 

Fig. [S] Inset shows the angular dependence of 



e„ (deg) 


FIG. 5. Angle dependent linewidth ratio 
The shadowed region shows the average with errorbar 
(1.50±0.02). Inset: Angular dependence of poAHi /2 for 
Pyo. 8 Cuo ,2 and Pyo. 8 Guo. 2 /Gu/Go at uj/2tt = 10 GHz. Solid 
lines are macrospin calculations. 


ATL^y^®^® (red) and AiL^l^ (blue) for Pyo. 8 Cuo. 2 . 
The data can be reproduced through a macrospin 
model 28, 2^ as shown in the solid curves, using similar 
magnetizations and isotropic dampings extracted from 
Fig. 2(a) and (c) (/rgMeff = 0.53 T, ai = 0.0114 
for the single layer, noMeS = 0.55 T, 03 = 0.0168 
for the trilayer). The inhomogeneous broadenings are 
negligible, shown in Fig. 2(c). For small enough On, 
the resonance field of the Co starts to fall below the 
expected macrospin value, as shown in the Supplemental 
Materials, Section cjl^. We take the angle at which 
this behavior appears (at 0 h 18°) to be the limit above 
which we have the greatest confidence in single-domain 
ordering of Co. 

In the main panel of Fig. [5] we replot the trilayer 
and single-layer linewidths for Pyo. 8 Cuo. 2 , shown in 
the inset, as the ratio . Because the 

inhomogeneous linewidths are negligible for the struc¬ 
tures (< 0.5 mT), the linewidth ratio for isotropic 
spin pumping would be approximated well through the 
ratio of the Gilbert damping for the two configurations, 
AH^y ^/= 1-1- Aosp/ai. We find that the 
linewidth ratio is in fact constant within experimental 
error, shown by the shaded region in Fig. [5] The 
blue dashed curve shows the expected behavior for 
anisotropic spin relaxation, assuming Af;. = 38 nm, with 
a marked decrease in the linewidth ratio for low angles 
0H. A best fit to these data returns Af;. < 1.1 nm. If 
we restrict our attention to field angles 0 h > 18°, above 
which we have conhdence in macrospin behavior of the 
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Co layer, the best fit is not changed greatly, with < 4 
nm, within experimental error of the transverse length 

Extrinsic effects, i.e. issues of sample quality, may 
play some role in the results. First, longitudinal spin 
relaxation lengths if equated with the spin diffusion 
length Asri . are inversely proportional to (defect-related) 
resistivity (SOj. However, four-point probe measure¬ 
ments of the resistivity of our Co (5 nm) films show 
25 /rf2-cm, comparable with the 18 ^Sl-cm reported 
in the room-temperature CPP-GMR experiment [Ilj, 
and therefore comparably long spin diffusion lengths 
should be expected. Second, we see that the spin mixing 
conductances of Pyi_a;Cua;/Cu/Co measured here 
are lower than those measured in Ref. [l^ , on structures 
deposited elsewhere. The most plausible source of the 
reduction is a more resistive Cu layer, which adds an 
additional resistive term ^ l^ | (2e^//i)fcuPCu to gig. 
Here however the bulk Cu properties should have little 
influence over either spin relaxation length and should 
not affect the anisotropy of spin relaxation strongly. 

Our estimate of A^^ in Co is consistent with a general 
observation that spin relaxation as measured in spin 
pumping/FMR is shorter-ranged than it is as measured 
in magnetotransport. In Pd and Pt, the characteristic 
relaxation lengths for dynamically pumped spin current 
are measured as 1-5 nm 15, ^ 


whereas in GMR 
they are closer to 10-20 nm[3^[37|. We suggest therefore 



that the quantities revealed by the two types of measure¬ 
ments may differ in some respect. For example, robust 
spin-pumping effects have been found in ferrimagnetic 
insulators such as yttrium iron garnet (YIG). These 
effects clearly have little to do with electronic transport 
in YIG, and their characteristic lengths would refer 
to scattering mechanisms distinct from those involved 
in GPP-GMR. A second possibility, alluded to in the 
review in Ref. 0, is that the room-temperature spin 
diffusion length of 38 nm in ll| is an overestimate due 
to technical issues of the GPP-GMR measurement in 
Co multilayers; the majority of such measurements in 
various ferromagnets show <10 nm[l^. Our results, in 
this scenario, may alternately imply that the short spin 
diffusion length observed in Py is not far away from that 
of Co. 

In summary, we have experimentally demonstrated 
that the spin relaxation in Co, as measured by non- 
collinear spin pumping, is largely isotropic. The 
estimated longitudinal spin relaxation length, < 2 nm, 
is an order of magnitude smaller than measured by 
magnetotransport but comparable to the transverse spin 
relaxation length. We acknowledge NSF-DMR-1411160 
for support. 
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A. Calculation of spin-pumping damping for noncollinearly magnetized, 
asymmetric trilayers 





FIG. 1. (a) Magnetization configuration of the asymmetric F 1 /N/F 2 trilayer, (b) An instant in 
which the spin polarization of is orthogonal to both mi and m 2 , is also orthogonal to 

mi and m 2 , (c) An instant in which the spin polarization of 7^“™^ is in the same plane of mi and 
m 2 . /.r^TV is also in the same plane of mi and m 2 

Consider an asymmetric ferromagnet / noble metal / ferromagnet (F 1 /N/F 2 ) spin-valve 
trilayer structure, shown in Fig. 1(a). The time-averaged magnetization of Fi is pictured 
along the film-normal, although it can take any angle with respect to the film-normal. We 
assume that Fi undergoes small-angle precession. F 2 is noncollinearly magnetized with 
respect to Fi, where 6 is the angle of noncollinearity or misalignment; 6 = 0 for parallel 
magnetizations mi = m 2 , where m, , z = 1 or 2, is the unit vector of the magnetization M* 
of Fj. The magnetizaton of F 2 is taken to be stationary. The spin current flows from the N 
spacer to each of the F layers Fi, F 2 arei^-: 

* ~n 

jN^Fi _ . mi) + ^mi x x mi (1) 
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* ~ti 

jN^F 2 ^ ■ m 2 ) + x x m 2 (2) 

where Hgj^f is the spin accumulation vector in the N layer, g* and gj'^ are the effective 
longitudinal spin conductance and transverse spin mixing conductance for Fj/N interface, 
respectively. Here the spin current vector denotes the direction of spin polarization, the 
direction of current flow always being normal to interfaces. The conservation of spin angular 
momentum, assuming spin-current conservation (negligible dissipation) N, gives: 

jN-i-Fi _|_ jN-fF2 _ jpump 


where is the pumped spin current from Fi into 


jpump ^ ^gl^mi X rill (4) 

Substituting Eq. o, o and dl]) into the continuity expression (|3]), we obtain a vector 
equation in terms of the vector spin accumulation To calculate the spin pumping 

damping enhancement, we seek solutions for in order to find the spin current flow into 
m 2 , which is absorbed by m 2 . 

The vector proportional to mi x riii, rotates in the plane with normal given by mi. 

Assuming a finite misalignment angle 6 between mi and m 2 , I^i™^ will oscillate between 
fully orthogonal to m 2 (Fig. lb) and canted away from orthogonality by 6 (Fig. Ic). We 
consider these two extreme cases during the precession of mi. In case 1 (Fig. lb), IP“™p is 
perpendicular to both mi and m 2 . In case 2 (Fig. Ic), IP“™p is in the same plane as mi 
and m 2 and has the largest longitudinal component along m 2 . 

In case 1, jN-i.F 2 gj,g gp parallel to In Eqs. ([T]) and ([2]) the first 

terms become zero and only the second terms remain. The solution of Eqs. ([I])-(I2D has a 
scalar form along the direction 5 orthogonal to both Fi and F 2 : 


jN^Fi{F2 } 


-tt 

3i{2} 


~tl I -tt 

^9i +92 


rpump 2 _ 
-'sl ) H'sN — 


4^ ( ^pump 

~tl I ~ti ' 

9i + 92 


(5) 


It has been shown previously that^'^ the dissipation of spin angular momentum due to a 
transverse spin mixing conductance g+ leads to an additional Gilbert damping term Aa = 
'-yfig^^^/AuMsd. With only dissipated, the spin-pumping damping enhancement can be 


expressed as: 


Afy = ■_ 

L-i.uLgp ± g2 jpump 


Aao _ . gf 

— Aao • 


~tl I -tt 
+ 92 


( 6 ) 
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with Aao = ■yhgl^/AnMsid. Eq. (E]) is identical to the collinear spin pumping case with an 
effective spin mixing conductance + (^ 2 ^)”^- 

In case 2, has only a component coplanar with mi and m 2 and /i^Ar)- In Eq. 
([I]) and ([2]) both terms need to be considered. The x and y components of Eq. ([2]) can be 
written as: 


dvr/” 



(-/^^Tv sin 0 + cos 9) cos 9 + sin 9) 

+fi'2/^sAr (7) 


The solution of Eq. ([7]) can be expressed as: 


/^sJV 

h-sAT 


47r/”(C'cos0-5sin0) 
AC-B^ 

47r/”(Asin0-5cos0) 

AC 


(8a) 

(8b) 


where 


A{9) =gl sin^ 9 -|- g\^ cos^ 9 -|- g\^ 

(9a) 

B{9) = (^1'*' — gl) sin 9 cos 9 

(9b) 

C{9) = gl cos^ 9 -I- ^1'*' sin^ 9 + g^ 

(9c) 


The spin torque is equal to the conponent of transverse to mi, or the component 

which is parallel to I^i“^. Thus the spin-pumping damping enhancement can be written in 
terms of the dehned misalignment-dependent quantities A{9), B{9), C{9) as: 


Aasp{9) 


jN^F^ . jpump 


=Aa, 


^pump y-pump 

g 2 {A sin^ 9 — Bsin 9 cos 9) + gl^{C cos^ 9 — Bsin 9 cos 9) 
^ AC -B^ 


( 10 ) 


It is easy to verify that at 0 = 0° Eq. ffTOj) recovers Eq. (|6]), same as the collinear spin 
pumping. 

Having treated the two special spin current orientations, we need to take the average of 
all the orientation possibilities. We refer to the calculation by Taniguchi, et ah-, that in a 
symmetric spin valve (Fi=F 2 ) the small-precession limit of averaged spin-pumping damping 
enhancement is equal to the arithmetic mean of damping enhancement with out-of-plane 
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(case 1) and in-plane 1^“™^ (case 2). The Eq. 13 in Ref.- can be simplified, at small 


precession angle, as: 


(i//2) sin^ e 


= A«o 1 - 2n 

1 — COS"' a 

which is the average of Aao and Aq;o[1 — psiv? 6 /(1 — cos^ 9)] (Eq. 5 in Ref.-). We apply 
it to the asymmetric spin valve condition: all the theoretical curves in the main text are 
calculated from the mean of Eq. ([6]) and Eq. OB. 

The theoretical curves in Fig. 4 and 5 of the main text are calculated using the routine, 
assuming = 38 nm for Co. The new estimation of (< 2 nm) in the manuscript takes 
the best value that hts the damping calculation to the experimental data. 


B. Values of g* and 


In this section we calculate the value of the two effective spin conductances. The trans¬ 
verse spin mixing conductance g^^ (Sharvin correction includedtserkovnyakRMP2005) of 
each interface can be calculated from the effective spin mixing conductance of Pyi_ 2 ;Cu 2 ;/Cu/Co 
structures and the comparison measurements of Pyi_a;Cu 3 ;/Cu/Pt and Co/Cu/Pt (Table [I]). 
For Pyi_ 3 ;Cua;/Cu/Co, the total spin mixing conductance can be expressed as: 


tt ~tt -tt 

^Pyi_:,Cu,/Cu/Co fi'Pyi_,Cu,/Cu S'Co/Cu 


( 12 ) 


For F/Cu/Pt (F=Pyi_ 3 ;Cua; or Co), the effective spin mixing conductance can be formulated 
as: 

1 1 1 
5'F/Cu/Pt S'f/Cu 5'Cu/Pt 

In the experiment the thicknesses of Pt are kept the same and we can treat as a 

constant. Solving Eq. flT^ and flT^ we obtain: 

11 11 


~tl tl tl tl 

5'Pyj_^Cu^/Cu 5'pyj_^Cux/Cu/Co 5'py^_^Cu:,/Cu/Pt 5'co/Cu/Pt 

11 11 

-+ 


~tl tl tl tl 

5'Co/Cu fi'Pyi.^^Cu^/Cu/Co fi'Pyi.^^Cu^/Cu/Pt 5'Co/Cu/Pt 


(14a) 

(14b) 


In Table UTI we list the calculated values of Cu^^/Cu ^ct/Cu- x = 0.1 and 0.3 

T I 

we take the linear interpolated values to evaluate gp^^ Cuj,/Cu/Pf addition, we also show 
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the values compensating the Sharvin correction, with l/g]^ = + l/2g^^ g^ = 15 

nm“^. 

Compared with previous measurementsl^iii, we hnd smaller values of g^^ /cu/co 
^Pyi Cu 2 ,/Cu a; = 0. However we argue that the spin mixing conductances of Co/Cu/Pt 
in Table [I and Co/Cu interfaces in Table UTl are reasonable, which ensures a good Co/Cu 
interface crucial for the study of spin relaxation anisotropy. It is also possible that a resistive 
Cu spacer contributes an additional resistance, (2e^//^)^cuPc^^^) to the right side of Eq. ([13]). 
To reduce the spin mixing conductance of Py/Cu/Co from 15.0 nm“^ in Ref.— to 7.6 nm“^ 
in Table (U one needs to take pcu = 16.8 fifl-cm. However we point out that this resistive 
scattering will contribute to both transverse and longitudinal spin conductance by the same 
amount, and the anisotropy of spin relaxation should not be affected. In practice, we take 
the effective interfacial spin mixing conductance into the model for the estimation of and 
and use the values of from Ref.—. 


(Unit: nm 

O 

II 

X = 0.1 X = 0.2 

II 

o 

CO 

II 

o 

tl 

^Pyi-r.Cu^/Cu/Co 

7.6 

5.6 7.3 

6.8 6.8 

tl 

5pyi_,Cu:,/Cu/Pt 

6.0 

5.0 

5.0 

tl 

^'Co/Cu/Pt 

7.9 




TABLE I. Experimental values of (effective) spin mixing conductance of Pyi_ 2 ,Cua,/Cu/Co, 
Pyi-^j-Cux/Cu/Pt and Co/Cu/Pt samples, extracted from spin-pumping linewidth enhancements. 


The effective longitudinal spin conductance g* can be expressed as^': 

1 _g^^ + g^ 1 

s' 2 s 1 St/ 4 S + g^jtunhitf/\i^) 


(15) 


In the hrst term, is the interfacial spin-up {spin-down} conductance. pg 

calculated by 1 /^RRl} = where is the electron interface resistance. 

We use the experimental value from GMR measurements: 2AR* = {AR^ + AR^)/2 = 1.04 
fH-m^ for Co/Cu^ and 1.0 fH-m^ for Py/Cu^. We can calculate that 2g'^g'^/{g^ + g^) = 26 
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(Unit: nm 

X = 0 

X = 0.1 

X = 0.2 

X = 0.3 

X = 0.4 

^Pyi_,CuUCu 

11.7 

8.6 

9.5 

9.1 

9.1 

-ti 

5'Co/Cu 

21.9 

16.2 

31.5 

27.2 

27.2 

tt 

^Pyi_.CuUCu 

8.4 

6.7 

7.2 

7.0 

7.0 

ti 

Sco/Cu 

12.7 

10.5 

15.4 

14.3 

14.3 


TABLE II. Sharvin-corrected spin mixing conductance of Py]^_j,Cua;/Cu and Co/Cu in- 

"T I 

terfaces, calculated from Eq. m- “^r”: interfacial spin mixing conductance compensating the 
Sharvin conductance of Cu layer, i = Pyi_a;Cux/Cu or Co/Cu. 


nm ^ for both interfaces. 

In the second term, Qsd has been expressed in Ref.- as: 

h 2(tV^ 


(16) 


where are the spin-up/down electron conductivity in F, h is the Planck constant and 
e is the electronic charge. Here we simply take = a'*' = cr/2 (a is the total electrical 
conductivity), which has also been done in Eq. (74) of Ref.-. Following this treatment, the 
term 2a'^+ a^) is replaced by a 12. Taking pco = 25 pfl-cm and ppy = 30 /ihl-cm from 
our four-point probe measurements and = 38 nm for Co^ and 4.3 nm for Py^ from the 
literatures, we calculate ( 7 srftanh(fi?/Af^) to be 0.18 nm“^ for Co and 8.3 nm“^ for Py when 
the F thickness is 5 nm; the large disagreement comes from the expected difference in Af'^. 
As a result, g* = 0.18 nm“^ for Co and 6.2 nm“^ for Py are obtained from Eq. (IT5l) and 
used to produce the theoretical curves in the manuscript. 

In the experiment, we do not hnd the anisotropic response of spin pumping predicted 
above. According to our model, the lack of anisotropic response can be explained best 
through a difference in the longitudinal spin conductance g* for Co/Cu, as this is the most 
dominant term in Eq. flTOll and sensitive to Af^. This is because in the experiments we choose 
the thickness of Co to be much less than 38 nm in order to examine the spin relaxation 
anisotropy. 

From Py to Pyi_a:Cux, we should expect that both g^j^d will increase due to a 

better conducting ability of Cu than Py. Af^ may also vary. However we emphasize that in 
Eq. ffTOj) . the anisotropy is dominated by g^ and gl'^ and not sensitive to g^. For example in 
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the angular dependence of linewidth ration for x = 0.2 (Fig. 5 of the main text), increasing 
fi'pyo sCuo 2 /Cu ^ factor of two will change the single-domain estimation of from 1.8±2.7 
to 2.1±2.8, still much smaller than the GMR measurements. Thus for simplicity we keep 
using the value of g* of Py for Pyi_ 3 .Cua; layers. 


C. Single-domain limit determination 



FIG. 2. Resonance peak of Co and Pyo. 8 Cuo .2 independently measured in Pyo. 8 Cuo, 2 /Cu/Co with 
Oh = 18° (a) and 9° (b). The resonance frequency of Pyo. 8 Cuo .2 are both 10 GHz. The resonance 
frequency of Co is adjusted so that the po-Rres of Co is equal to Pyo. 8 Cuo. 2 - Dashed curves show 
the theoretical prediction of Co resonance signals. 

To determine whether the Co layer is in a single-domain state in Pyo. 8 Cuo. 2 /Cu/Co when 
the Pyo. 8 Cuo .2 layer is at resonance, we have measured the FMR signal of Co at different Oh- 
First we measure the FMR signal of Pyo. 8 Cuo .2 at one angle and determine the resonance 
held po-^res- Next we adjust the frequency so that the Co FMR signal can be measured at 
the same held. Then we compare the lineshape with the macrospin model prediction!^. In 
Fig. 2(a), when Oh = 18° the resonance held po-f^res for Pyo. 8 Cuo .2 is 0.53 T at 10 CHz. 
For Co, the resonance held is located at 0.53 T for 14.8 CHz. The macrospin model for 
angle-dependent FMR shows u:/2'k = 14.8 CHz, identical to the experiment, showing the 
Co can indeed be treated as a macrospin. However for Oh = 9° (Fig. 2b), we hnd that the 
Co resonance is located at 15.0 CHz, quite diherent from the macrospin prediction of 12.3 
CHz. To see the diherence more clearly, we plot (dashed lines) the macrospin prediction for 
both Pyo. 8 Crio .2 and Co resonances at 15.0 CHz, based on the magnetizations and linewidths 








measured from perpendicular FMR. The Pyo. 8 Cuo .2 peak matches with experiment. The 
calculated Co peak deviates from experiment in both resonance held and linewidth. Thus 
we determine the single-domain limit of Oh to be somewhere between 9° and 18° in the 
sample. The upper bound 18° is used in the manuscript for the single-domain limit. 
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